function ordertest()
  GRID_SIZE = 64;
  L = 0.16;
  Z_SIZE = [32 64 128 256];
  tfinal = 20;
  TIME_STEPS = 20;
  GAUSS_AMP = [0 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.25 0.30 0.35 0.40 0.45 0.5];
  GAUSS_SIZE = 0.1;

  res = zeros(length(Z_SIZE), length(GAUSS_AMP));

  for z=1:length(Z_SIZE)
    for g=1:length(GAUSS_AMP)
      disp(sprintf('*** z-size: %d,  amp: %.2f', Z_SIZE(z), GAUSS_AMP(g)));
      r = h2amz(GRID_SIZE, L, Z_SIZE(z), TIME_STEP, TIME_STEP * TIME_STEPS, GAUSS_AMP(g), GAUSS_SIZE);
      r(TIME_STEPS)
      res(z, g) = r(TIME_STEPS);
    end
  end
  save('ordertest.dat', 'res');






end
